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1. Introduction 


The differential equations studied earlier often possess solutions 
expressible in terms of elementary functions such as polynomials, 
exponentials and trigonometric. However, many equations arise from 
physical problem do not have solutions that can be so conveniently 
expressed. 


In this chapter, our goal is to obtain representations for ODE solutions 
as power series. One of the major benefits of using power series Is 
that it will allow us to solve many ODE involving variable (non- 
constant) coefficients. 


2. Review of Power Series 


Power Series 
The power series (in powers of x — Xo) 1S an infinite series of the 
form 
on 
im 2 
(1)... Ya, (x-%) =a, +a, (x -X,))+a,(x- X))o te, 
m=() 


A),4\,Qy,... are constants, called the coefficients of the series. xo is a 


constant, called the center of the series, and x is a variable. 
[f in particular xp= 0, we obtain a power series in powers of x 


| a 2 
(2)... ya,x =A) t+QX+a,x” +e 
m=() 
assuming that all variables and constants are real. 


Remembering how Taylor series work will be a very convenient way 
to get comfortable with power series before we start looking at 
differential equations. 


Taylor Series 


If f(x) 1s an infinitely differential function then Taylor series of f(x) 
about x=Xp 1S: 


00 (m) 
foay= yA S05," 


Recall that: 


FO (x) = F(x) & f(x) =m" derivative of f (x) 





= 


UCSI Example 1: 


DY tvs Taylor series for f(x)=cos(x) 


According to Taylor series f(x)=cos(x) can be expressed: 

£ (x) =e08(x) 
f® (x) =—sin(x) 
| f°) (x)=—c0s(x) 
_ f°'(0) n in (0) ff ‘i (0) x2 4 f° (0) eee | i (x) =sin (x) 


= Xr 
f (x) =cos(x) 


© (mm) 
cos( x) = 2d Fy —X,)" 








0! l! ‘an 3! 


oo 4 6 8 
=! 49-2 40440-40424... 
0! 2! A! 6! 8! 





f (0)=1 
f® (0) _¢ 
f (0) | 
f© (0) =i 
ft (0) _| 


Example 2: 
Determine Taylor series for f(x)=e* about x,)=0. 


Solution 
This 1s probably one of the easiest functions to find the Taylor series for. We just need to recall 
that. 


f” (x ) oe Bs 0.1.2... 
and so we get, 
f” (()) =| i= (). l. i, 7 


The Taylor series for this example 1s then. : 
e =? 


Familiar examples of power series are MacLaurin series. 





3. Theory of Power Series Method 


Convergence and divergence 
CO , ; 
bo » 4, (X— Xp) SH, FOAL PAAR AKA) Pee, 
m=0 


If the power series equals a finite real constant for the 
given x, then the series 1s said to converge at x (the 
series 1s convergent for the given x). 


If the series does not converge at x, 1t 1s said to 
diverge at x (it is divergent for the given x). 


Note that polynomials are also power series. 


Convergence Interval & Radius of convergence 
Theorem 


= > iy (X—%y)" = Ay i e—-h Teese. 


m=(0 


be a power series ,then there exists some number O<Rsece 
,called its radius of convergence such that the series is 


convergent for _, (x=) < R 
and divergent for — (x = “> >K 


The values of x, for which the series converges, form an interval 
called the convergence interval. 


Every power series has an interval of convergence. 


Radius of convergence 


What is R? 


The number R is called the radius of convergence. It can 
be obtained as, 


R= l 
OR lim ¥/la,, 








Radius of convergence 


Example 1 The useless Case 1 of convergence only at the center 
In the case of the series 


oo 
> mix™ =14+ x4 2x7 + 6x ++, 








m=0 
: 1S mt] 
we have a,, =m! andin R=1/ lim |——|, 
mc] gq 
it 
a m+1)! 
@m+1 _ (im +I)! yt em+1—>0 as m1 —> OO. 
a, m! 


Thus this series converges only at the center x = 0. Such a series is useless. 


Example 2 Case 2. The geometric series 
For the geometric series we have 


| _ ‘ ae 2 
rag ~-1Lt+x+x° 4+: (jx{<1). 
m=0 
In fact, a, =1 for all m, and aay tim “m) we obtain R = |, that is, the 
HEP OO ay 








geometric series converges and represents 1/(1 — x) when | x!<1. 


Example 3 Case 3 
In the case of the series 


ora ” x? 
ea) =l¢x+O+-- 
| 2! 











m=() 
, fa 
we have a_ =1/m!. Hence in R=1/ lim|-“*! 
. mr] Qh 
a I/(m +1)! | 
mat _ 1/( oil oy aS M1 —-> 0, 
a, 1/ni! m+] 


so that the series converges for all x. 


Exercise 1: 


Determine the radius of convergence, R, for the 
following power series: 


“mx” 
» ee 


m=0 





(-1)"x" 
aro 


4. Existence of Power Series Solutions; 
Real Analytic Functions 
Definition: Ordinary and Singular Points 


y''+ p(x) y'+q(x)y = 0 (3) 


A point x, 1s said to be an ordinary point of the DE 
(3) if both p(x) and g(x) are analytic at xy. A point that 
is not an ordinary point 1s said to be a singular point 


of the equation. 


What is a real analytic function ? 
A real function y(x) 1s called analytic at an ordinary 


point x = x, if it can be represented by a power series 
y P (x ~ x, )" with radius of convergence R > 0. 
m=(0 


. . . Examples 
Ordinary and Singular Points 


(a) The differential equation xy" + (sin x)y = 0. A point at x = O 1s 
ordinary point since it can be shown that g(x) = (sin x)/x possesses 
the power series form: 


+ + 
3! 0ST! 


that converges for all finite values of x. (R=) 


(b) The differential equation y” + (in x)y = O has a singular point at 
x = O because g(x) = In x possesses no power Series 1n x. 
(Note: /n(Q)=-00) 


(c) The singular points of the equation (x - J)y" + 2xy'+ 6y =0 are 
the solutions of x* - 7 = 0 or x = 4. D(x)= 2x/(x° - 1)=00 and also 
q(x)= 6/(x° - 1)=0 for x=+/. All other finite values of x are 
ordinary points. 


(d) xy" - 2xy' + 2y = 0 and xy" - 2y = O has a singular point at 
x = O since p(x)= 2/x=oo and also g(x)= 2/x*=00 for x=0. 


4. Existence of Power Series Solutions: 
Real Analytic Functions 
Existence of the power series solutions 


y"'+ p(x)y'+g(x)y =0 = 
The power series solution to DE (3) will exist and be 
unique provided that p(x) and g(x) are analytic in the 


domain of interest (interval). 


Thus, 1f p(x) and g(x) have power series representations, 
then (3) has power series solutions. 


5. Power Series Solutions of ODE 
Idea of the Power Series Method 


For a given equation 
yit p(x)y'tq(x)y = 0 “ne 
where p(x) and g(x) are polynomials. We assume a solution in the form of 
a power series with unknown coefficients, 


oo 
yn dq” =) +a,Xx+ ax fees ...(4) 
m=0 


5. Power Series Solutions to ODE 


By term-wise differentiation we obtain y’ and y" 


CO 


y= > 4," =A) FAX+ A,X FOX FH ecccccees. (4) 
m=0 
SO, y= >» ma,,x" =A, + 20,X+4+30,X" trees. (4a) 
m=| 
= » m(m—l)a, x" ~ = 2a, +3.2a,X+4.304X° veces. (4b) 


m=2 


5. Power Series Solutions to ODE 


Step of the solution 


> Substitute equations (4), (4.a) and/or (4.b) into 
original DE. 


> Collect like powers of x. 


> Equate the coefficients of each occurring power of x 
to zero. 


> Hence the equation can expressed in power series 
form. 


Example 3: 
Solve p'- y=0. 
Solution y= Ye, x" =a, +a,xta, x2 +a, 4. (3) 
y= S'ma,x”” =a, + 2a,x+3a,x* ++ (3a) 


ni=!l 


Insert these into the given differential equation 
(a, + 2a,x +3a,x° +---)— (ay tax t+ A)x° ++)=0. 
Then we collect like powers of x, finding 
(a, ~ ay) + (2a, — a,)x + 3a; — a,)x° +--+ =0. 
Equating the coeffictent of each power of x to zero, 
a, — Ay =0, 2a, — a, =0, 3a; —a, =0,:: 
“~ 9 =-2.% .., 


GQ, -da q,=—--= 
‘ii. = = 
| (jo **? 3 - 
2 2! 3! 


Example 3 cont'd... 


a a a, a 
a, =@,a,=~7=—, 0, = SHS ee 
2 2! 3 3! 
With these coefficients, y= a, + a,x +a)x° +.a,x* +--- becomes 


a, Ay 3 
Y=A) tayx+—x? +2 x? +. 


3 
x x? 
malloc 44 


—+-:+|  @a@m Power series 
3! tt 


= Aye. 


We can conform the result by using the basic technique: separation 
of variables > integration. 


Example 4: 
Solve p°+yp=0. 


oo 
Solution y= > a,x" =ag + axtanx* +a3x° +--, (3) 
m=0 
- 
y=) max”, (3a) 


m=! 


oo 
7 = S mm — Da x”? =2a,+3-2a,x+4- 3a,x* +- (3b) 
ni=2 
By inserting these into the equation we obtain 


(2a, +3-2a,x+4-3a,x* poe) (0, + a,x +a,Xx* +---)=Q. 
Collecting like powers of x 
(2a, +a.) +(3-2a, +a, )x+(4-3a, +a,)x° ++ =0, 


Equating the coefficient of each power of x to zero 
2a, + ay = 0 coefficient of x” 


3- 2a, +a = Q coefficient of x’ 
4 ‘3a, +a, = O coefficient of x’. 


Solving these equations, 


Example 4 cont’d... 


dg and a, are arbitrary. 
With these coefficients, y =a, + a,x +a,x° +a;x° +++» becomes 


a a ad a 
Y=) tax, —-2 x? — x? + 2 xt +t? +. 
2 3! 4! 5! 


—qQ pee +a a 
OO or Fy 31 5 
NS SS >—— 


y = dy.coS(x)+ a;.sin(x) 


Note: 


~~~ 1. It Is not always necessary to get the power series solution in 
the form of familiar functions! (trigonometric, exponential, 
etc), 


2. In example 1 and 2, the series solutions are done in open 
form (the series are expanded). Later we will see how the 
solutions are done in closed form. 





Closed form 


Exercise 2: 


Solve: y’+2xy =0 


By using power series method! 


6. Operations on Power series 


(1) Termwise Differentiation 


lf Y= > nx” is convergent, then 
m=0 


y = > ma,.x"" and y= S" m(m—la,.x"” 
are sie convergent. - 
(2) Termwise Addition 
If Yay" and Db are convergent 


in the same domain x, then the sum__ > (a,, +8,,)x” 
also converges in that domain. “ 


(3) Termwise Multiplication 


Two power series may be multiplied term by term. 
If >. A,X” and Ds Pn® are convergent. 
m=0 os 


Then the series obtained by multiplying each term of the 
first series by each term of the second series and collecting 
like powers of x is also convergent. 
» (a,b, +ab,,+...+a,b,)x" 
m=0 
(4) Shifting Summation Indices 
This is technically best explained by example (we will see 
in the last example). 


Examples 5 


Solve by power Series: y"- 2xy = 0. 


Solution: 
x = 0 1s an ordinary point of the equation. From the above Theorem, 


we can find the solutions of the power series form: 
y= Dex” 
m=0 


— a >» me,,x" 


m=l 


y= » m(m—l1)c,.x"~ 
m=2 


y'-2xy = > m(m —l)c,x"* - ae ie = 0 
m=2 m=0 


Therefore, 


Examples 5 cont'd... 


m=2 


y"—2xy = > mm —lc,x"* - » 2¢,,%"" = 0) 
m=0 


Now we will need to shift the first series down by 2 and the second 
series up by | to get both of the series in terms of x*. Let k = m - 2 in 
the first series and k = m + J in the second, we have 


> (k+2)k + Vep..x" — ) 2c, x" =0 
k=0 k=] 


Next we need to get the two series starting at the same value of index k. 
The only way to do that for this problem 1s to strip out the k=0 term. 


(2)(1)c, (1) + S (k + 2)(k +1)c,,,x* - S2¢,_,x! =0 


2c, + S[k + 2)(k +1)c,,, —2c,_, * =0 


k=] 


Examples 5 cont'd... 
We must then have 
20=0 Do=0 for k=0 
(k+ 2)(k + [)cp4.-2c,.; = 0 fork = /, 2, 3, ... 





and 
, _ 2Cp_] 
k+2 (k+2)(k +1)? fork = J, ve a 
This leads to recurrence relation: 
7 
4 
7.6" 7-6-4-3 
2Cy 9) 
So a = 0) 
_ to 
C 
C= 73 2 23 
ii (9 = 7 5 =O ope £0 
2c 9.8 9-8-6-5-3°-2 
at 
5-4 2 3 
CQ =—c = ———_¢ 
ae 1010-9. > 10-9-7-6-4-3 | 
6 3 0 
6-5 6-5-3-2 


Examples 5 cont'd... 


The solution is now: 
CO 
Y= DCX” 
m=0 


Vy b0ieteod” beak +egn $eex” Fe 
2 52 


6 7 
———— Cox 
One oaZ 


7 vs se 4 
=Co +c;x+04+—cox +—c]x 1x 
3-2 


+————_¢ 
7-6-4-3 


2 a y) ae 
=v fo ee +Cy it 
2 6-5-3-2 44 7-6-4-3 


y =cyy,(x)+¢,y,(x) " 


Examples 6 


We will repeat Example 4, but now the power series will be 
expressed in Closed form. 


Solve it by power series: y"'+y=0. 
CO 
y= a,x" 
m=0 
i > mc,,x = >» ma,,.X" 
m=0 m=1 


= > mm- la, x" ~* = > mm- l)a, x" 
m=0 m=2 
Substitute them into the equation, therefore: 


y m(m—l)a, x" + 3 as = 
m=2 


Examples 6 cont'd 


> m(m—1)a,,x""~ + ys a,x" =0 
m=0 


m=2 


— 


2D (k+2)(k+1)ay..x° + > 4x = 0 
k=0 7 


OO 


S[(k + 2)(k +Da,,> +a, ]x* =0 


k=0 
(kK+2)(k +1)a,,, +a, =90 For k=0,1,2....... 

Since the problem is 2" order DE, we a, 

will let a, and a, be unknowns O42 ~ (k+2\(k +1) 


(arbitrary), thus: 


Examples 6 cont’d 


Ay 


3 For k=0,1,2....... 
’ (k+2)(k +1) 


With these coefficients, y = ay) + a,x +a)x° +a;x° +++ becomes 


a a ad d 
Y=, tax, —~—o x? ~ bx? 4 Ext +t? 4+... 
2 3! 4! 5! 


x? x4 - x 
y=a, l-—+-—-— +: ]|+a] x-—+—— + 
. 21 4! 3 OCS 
| if 
yY = dy.COS(x)+ a).sin(x) 


Exercise 3: 


Solve this ODE using power series: y" + xy’ + y =O 





Thank you 


To be continue ... 


